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Abstract 

The nonrelativistic guiding center motion of a charged particle in a 
static magnetic field is derived using the Hamiltonian formalism. By 
repeated application of first-order canonical perturbation theory the 
first two adiabatic invariants and their averaged Hamiltonians are derived, 
including the first-order correction terms. Other features of guiding center 
theory are also obtained, including lowest order drifts and the flux invariant « 
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CANONICAL PERTURBATION THEORY 


This is an attempt to formulate and solve the nonrelativistic guiding center 
motion of a charged particle in Hamiltonian terms. In other words, since 
guiding center motion may be viewed as a perturbed periodic motion, we try 
to apply to it methods of celestial mechanics (suitably modified) which are 
designed for motions of this type. 

One such method (often associated with the names of Poincare and Von 
Zeipel^ is the following. One begins by obtaining the zero-order 

Hamiltonian — that is, the limiting Hamiltonian for the case of the 

perturbation tending to zero — and solves the associated Hamilton— Jacobi 
equation for action-angle variables (j, fl ). This should be possible since 
it is given that the £ — » 0 limit is both periodic and soluble. 

The Hamilton-Jacobi equation gives the generating function W of a 
canonical transformation to new variables (p^> 0^) that include the action- 
angle pair. In the system described by alone, all new variables are 

constants of the motion, except for fl which is linear in time and increases 
by unity each period. This transformation is now applied to the full Hamil- 
tonian (with finite £ ), so that the new variables vary slowly , except for 

XL which is now approximately linear in time. 

The prescription next calls for a near-identity canonical transformation 
(P ± . %) * ( 


with transformed action-angle variables (J*, J~L*) and with 
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H*( p , £j ) = H* (0) ( f> , <| ) + £h* ( 1) + ... (l) 

the new Hamiltonian, such that H* does not contain f\* . If this is 
accomplished then J* is a constant of the motion, termed adiabatic invariant 
in guid ing center motion and sometimes called a third integral in celestial 
mechanics (the terms are not synonymous and their difference will presently 
be noted). Further solution of the problem by means of H* does not have to 
contend with (J*, CL*) and thus involves two variables less than the origi- 
nal problem, corresponding to a reduction by one dimension. 

Let the near-identity transformation have a generating function 

S(p,q) = Lpili ^ JZ SW( P- ' I > (2) 

k*l 

Then if the time t is not explicitely involved, the old and new Hamiltonians 
are equal, giving 

22E k B w (p.*Le."'»s w /^.a) = 

k=0 

- JZ e k H* (k) ( p . 5 + £ e" 'asW/'sp. ) 

k=0 

Expansion gives, for 0( iP) 

H (0) ( h . q ) = H* (0) ( p , q ) (3) 

which is acceptable, as H ^ does not contain ii , the conjugate J 
being a constant of the motion. For the terms of order £ one gets 
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H tx) + 22 ( ^ (l)/ ^i )( ^ E(0)/ V “ 

= H^ 1 * + w 

If this is not to be a partial differential equation for , only 

one term containing may be allowed . In celestial mechanics, where 

the perturbation is "small" and is expressed by terms of various orders 
in £, added to this is accomplished by allowing ^ to contain 

only J : 

jj(o) _ jo)/ 2TI ( 00 = const,.) (5) 

Then ( 4 ) gives 

(c0/2h )(DS^/^CL ) + = H*^ ( 6 ) 

If (as is assumed) any dependence on fL is via periodic functions, then 
any function F containing jTL may be resolved into an averaged part 

1 

< p > = 5 p dn - ^ 

0 

and a "purely periodic" or "oscillating" part, with zero average, 

(P) = P - <F> 

' osc 

Now H*^ has no dependence on jfL while (Qs (l V^ ) is purely 
periodic — any part of independent of ft is eliminated by the 

differentiation — so 
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H*^ = <( (9) 

s (l) = - (2JT/0J) j"(H (1) ) osc da (10) 

which gives the transformation to first order. 

In guiding center motion the situation differs somewhat, since the 

perturbation is not necessarily small, but rather is what is known as 

slow or adiabatic (this leads to the difference between third 

integrals and adiabatic invariants). In an adiabatically perturbed system 

some of the variables (collectively denoted as q and p ) are "slow" and 

have the property that for any function F appearing in the calculation and 

depending on them -.(but not containing than as explicit factors), 

or F/'d P a are of" the order of £F . Slow explicit dependence on time 

is also possible, but will not be discussed here; it has been investigated by 

Gardner^, in an article discussed in more detail by Contopoulos^. Gardner's 

method differs somewhat from the one described here and further developed by 

( 6 ) 

Stern t in that it uses a separate canonical transformation for each order 
in t o 

In adiabatically perturbed motion may depend on variables other 

than J , provided that they belong to canonical pairs of which at least one 
member is "slow". Corresponding terms in (4) are then shifted to higher orders 
and one again ends up with (6), (9) and (10). Higher-order results differ 
somewhat from those derived for "small" perturbations, but in this work only 
the first-order correction terms will be considered . 
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Another complicating factor that only arises with orders higher than the 

(k) 

first is caused by the fact that the functions S' representing the 
transformation depend on "mixed" variables — old coordinates and new 
momenta. To express new variables in terms of old ones, or vice versa, 
further untangling is needed; this may be avoided by using the direct form 
of near-identity canonical transformations^^ (10)^ Here such methods are 
not required, however, for in calculations correct to the first order 
it is always permissible — due to the transformation being a near-identity 
one — to replace new variables by old ones in first-order correction terms. 

For instance, from the transformation relation 

Pi = = p i + 

and from (6) and ( 9 ), the adiabatic invariant will be 

j* = j - e W l} A)n + o( e 2 ) 

= J + £ (2^T/C0) ( q] ) osc + °^ 2 ) 

= J + £ (2JT/Cd) ( [p, q] ) qsc + 0( £ 2 ) (ll) 

Because the first-order correction of J* is purely periodic, the 

2 

long-term average of J is conserved to order £ .It may be shown from 
this^ that the total variation of J over a time period of order £ 1 is 
of order £ , leading to the well-known property of lowest-order adiabatic 

invariants, namely, that in a system undergoing a finite perturbation their 
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total variation nay be made arbitrarily small by stretching out the pertur- 
bation over a long enough time. 

The arguments of H* evolve on a time scale one order slower than 
that of the zero-order gyration of H . If however the motion represented 
by H* is also periodic in its lowest order, the preceding procedure may be 
repeated, leading to additional adiabatic invariants (this was also pointed 
out by Gardner). A particle in the earth's radiation belt, for instance, may 
have three independent periodicities, each associated with an adiabatic 
invariant, and these invariants should be derivable by repeated application 
of the preceding routine. 

Almost all perturbation methods contain a "smallness parameter" 

£ ^ 1 , but in adiabatic perturbation theory & is somewhat artificial^. 

It may be introduced by rewriting 0(£ ) terms (e.g. for some "slow" 
variable q g ^) as 

QF/0q si > £ (cF/ c ( tq 3 .) ) (12) 

but (since no definite numerical value of £ is available) it must be removed 
by the reverse process before the final result is obtained. The only time £ 
is used is in grouping terms according to the order to which they belong. In 
what follows there will generally be no question about the order to which a 
term belongs, and therefore we may avoid any use of £- altogether , although 
first-order terms will sometimes be labeled by superscript "(l)" . 




- 7 - 


PFF.T.TMINARY TRANSFORMATION OF THE HAMILTONIAN 

The guiding center motion of a charged particle in a magnetic field is 
strongly related to the structure of field lines and one would expect this 
structure to enter somehow into the calculations This is done 

via the Euler potentials ^ 11 \ oC , ^ ) > related to the field vector B by 

B = V<* * V^> 

Each field line then may be viewed as the intersection of two surfaces 


= constant 



constant 


f 


0 

0 


.A 

and^is consequently labeled by the parameters ( °^q* j o ° Associated with 
( , p> ) is a vector potential orthogonal to B , which will be 

henceforth adopted „ 


Euler potentials may be introduced into the dynamical problem of charged 
particle motion by choosing ( ©<• , f ) as curvilinear coordinates of position. 
For the third coordinate one may choose s , the distance measured along field 
lines from some given reference surface. In curl-free magnetic fields the scalar 
magnetic potential If would probably be a better choice, but this was not 
investigated. 

Gardner^ has carried out the first steps in transforming the Hamil- 
tonian and we shall start with his transformation, continuing from there 
according to the prescription outlined earlier. Our choice of subscripts 
follows Gardner's, but his use of £ is not adoptedo 
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Let the vector x give the cartesian coordinates of a particle of mass 
m and charge e , and let TT be the canonical momentum vector conjugate 
to x . The non-relativistic Hamiltonian for the particle's motion in a time- 
independent magnetic field is then 

n - (e/c) c< Vjl } 2 (14) 

A canonical transformation to new variables (P, Q) , generated by the function 

F(P, x) = o^(x) + P-l f (x) + P 2 s(x) (15) 

transforms to a new system in which the canonical coordinates equal 

(<* f ^,s).It does not, however, separate the rapid gyration around 

field lines from other motions, is achieved by a canonical transformation 

(7) 

proposed by Gardner , generated by 

F = P^oC + P i f + P 2 8 “ (c/e) P 5 P x (16) 

This gives 

H = (l/ 2 m) | P 3 Vo< - (e/c) + P 2 Vsj 2 (17) 

and 

= (c/e) P 1 + 0^ 

= + (c/e) P 3 

3 " S 

l 

In a homogeneous field with B = £ B ( £ unit vector ) this is readily 
solved ^ 11 ^ by choosing 
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oC = X ? =By S = z (19) 

leading to 

H = (l/2m) [p^ 2 + (e/c^B 2 ^ 2 + p/ ] (20) 

The variables (P^, Q^) appear here in the same form as the canonical 
variables in the Hamiltonian of the simple harmonic oscillator; therefore > 
they vary periodically and represent the gyration. On the other hand, by (18), 
c P^/e and may be identified as the (constant) Buler potentials of the 
guiding field line, while P 2 represents the momentum component parallel to 
the field and is also constant. 


Unfortunately, this solution cannot be immediately extended to the 

general nearly-homogeneous field, because the vectors ( V* , Vf , Vs ) are 

usually not orthogonal. The non-orthogonality of ( ) can be remedied 

2 

by adding to the generating function P a term ?>P^ , where 

/ 

> = (c/e) ( Vo<*v¥>)/ ( vp ) 2 (21) 


The non-orthogonality of Vs to the two other vectors is harder to correct. 
One might hope that some scalar s' replacing s would have the required 
property, but this is not so, except for special cases (e.g. curl-free fields, 
where the scalar potential may serve as s' ). For if s' existed such 
that Vs'-V°C and Vs'- Vp* Both vanished, then a scalar K must exist 


such that 


implying 


B = 


r 


Vs' 


B • ( V x B ) = 0 


which is not always satisfied. 
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However, if the field varies slowly in space, so that quantities such as 
VV<^ and are 0(£ ) , it is possible to define a scalar <T such 

that V(T is orthogonal to V<X and within 0(£ ) correction terms. We 
take 

(T = s- o(a-^b (22) 

where 

a = B" 2 B • ( Vs X 7&) 

(23) 

b = B" 2 B • ( X 7b ) 

giving ( B = B / B ) 

B t* Vs - aV«^ - bV^ (24) 

and hence 

V(T = S - o( Va - jWb 

- B + o(e) (25) 

In this calculation the field B ' and associated quantities such as Vo< , 
Vji , a, b and } are all 0(l) , while quantities derived from them by 
differentiation — V Vo(, Va , etc. — are all 0(€. ). However, o( 
and ^ appearing explicitely in eq. (25) are not of order £ ^ but 
also 0(l). The difference between them and 0(l) field quantities such as 
B is that field quantities depend slowly (or adiabatically) on position vari- 
ables, whereas o( and P depend on such variables in a "regular” (non-slow) 
fashion. This is evident from (18) for the variables generated by (16) and 
also holds in what follows. 
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To accomodate the non-orthogonality of the basic vectors we let the 
initial canonical transformation be generated by 

F(P, x) = o<P 3 + pP x + <T P 2 - (c/e)P 3 P 1 + i^P 3 2 (26) 

Application of the transformation equations gives 


0< St 

(c/e) P x + 

(Qj+ >P 3 ) 

(27) 

t - 

*1 + 

(c/e) P 3 

(28) 

cr » 

Qo 


(29) 


Also, the canonical momentum vector is (gradient taken in z space) 

7 T = vf = p 3 vo 4 + PjVf + p 2 v<r - iP 3 2 vA 

Substituting everything in (14) then yields 

H = (l/2m){p 3 V<< - (•/c)(Q 3 + >P 3 )7^ + P 2 V<T - ± | 2 (30) 

The vector being squared is the (non— canonical) momentum m v of the particle 
and the first two terms in it may be regarded as the zero-order component 
m v^ of this vector orthogonal to the magnetic field. If one defines a 
vector T orthogonal to both B and 

T = VoC (e/c)>V^ (31) 

then this part may be written in terms of orthogonal vectors as 

m v[ 0) = P 3 T - (e/c) vf (32) 
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THE GUIDING CENTER EXPANSION 
The magnitude of T is 

|T| = B/|vf| (33) 

which with (30) and (32) gives 

H = (l/2m) £p 5 2 B 2 /(7f) 2 + (e/c) 2 Q^Vf ) 2 + P 2 2 (7(r) 2 | + 0(£ ) 

(34) 

This resembles (20), but there exists an important difference: the variables 
(P^, Q^) appear not only as explicit factors but sIbo by (27)-(29), in the 
variables ( o( , p f ) giving the locations at which B and are 

to be evaluated. To eliminate (P^, Q^) from these variables one must 
expand any function of position around the "guiding center" point, denoted 
by subscript "c" and defined by 

( °<o' Pc’ <r 0 ) - (<*]/•.«!. «2 ) (35) 

By (27)-(30) , the guiding center expansion for any function f of position 
may be written 

f(o^,p. <r ) = ( exp D ) f( oi Q , f c , <T) 

where D is the operator 

D . (« j + Ap 3 )('»/'»<* c ) + (oP J /.)('i)/'»p 0 ) (36) 

and where (exp D) is defined by the Taylor series of the exponential 
function. In what follows, all functions depending on position will be 
assumed to be evaluated at the guiding center, unless otherwise is stated. 
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To understand the significance of D , let us define (in guiding center 
variables, subscript omitted) the operator 

(37) 

In addition, let the angular frequency of gyration be defined as 


a) = B e / m c 
and the gyration radius as the vector 



(38) 


(39) 


with v^ defined as in (32) but evaluated at the guiding center. One 
then finds (it is best not to simplify the cross products below) 


(f -vj 


- (“ B)' 1 



. \ (B X VoOC’D/'DcO 


+ 


+ 


(Bxvp/Df)j = 


(40) 


Two comments may be made here. First, in order for the expansion to be 
valid we require 

D « 1 


i.e* 

( ? *VJ«1 


(41) 


In other words, 2 is to be much smaller than the scale on which the field 
variables change. This is known as Alfven's criterion for the validity of 
the guiding center approximation and provides a more precise formulation of 
the statement made earlier that "the field varies slowly in space." 
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Secondly, the difference between 2‘ V A and P* V is small, since B 
and VO" are almost parallel; it is given by 

( 2 - vg ) 'V'acr = to -1 (v<r x b) • 'D/'fccr 

» o(e) 'a/'zxr (42) 

In the present calculation we only retain terms of order £ and therefore 
in such terms the difference between 2 m V± and «P*V may be ignored. 

Application of the expansion to the Hamiltonian (30) then gives 

2 

H = (l/2i) [mv^+ P 2 V<T + mv^} + 0(S 2 ) (43) 

where 

m - m j>»Vv^ + *2$' VV<r “ = 0(e) (44) 

and where all variables are defined at the guiding center. 

ZERO ORDER ACTION - ANGLE VARIABLES 


We now apply to the system the transformation derived from the Hamilton- 

Jacobi equation in the limit £, — ► 0 . In this limit both B and ( V^) 

are constant and the Hamilton-Jacobi equation derived from (34) separates 

into two equations, representing motion parallel and perpendicular to the 

field. Of these the former gives constant motion along the field's direction, 

while the latter resembles the Hamilton-Jacobi equation of a harmonic oscil- 

lator v A o Of course, when we actually apply this transformation to the 
2 2 

given system, B and {V $ ) must be allowed to vary slowly again. 
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If the new variables are denoted by (p, q) , with action-angle variables 


(Py 9^) ® (45) 

then the generating function of the transformation is 


w(p, a) 



(j$ e/jt c) 


(eS/c) 2 Qj 2 J 


£ 


+ 


PlS. + ^2 ^2 


(46) 


where 

5 = (Vf) 2 / B 


(47) 


and where all variables are evaluated at the "mixed coordinate guiding 
center" r' , defined by 

( * Vo * " ( c P]_ /• * \ > ^2 ^ t 48 ) 

We then have 

= < 7)W/'7>Q 5 » [ (JS e/fl c) - (e5/c) 2 Q^ 2 J (49) 

from which the first part of (34) is 

(l/2m) [ P 3 2 B 2 /(Vf ) 2 + (e/c) 2 Qj 2 (Vf ) 2 ] * W2JT (50) 


Strictly speaking, this may not be directly substituted in the Hamiltonian 
since (*) and other quantities are here defined in mixed variables. However, 
we have for such quantities 



- 16 - 


CO (r' c ) = CO(r c ) + (c/e) - P x ) ( / D<«>/Oc< ) + ... 

- <^(e c ) + 0(e 2 ) (51) 

since to/'Soi , (P^ - p^) and similar quantities are all 0( 8 ) . Because 

the present calculation ignores second-order terms, the difference between 
"mixed" and "true" guiding center variables (either "new" or "old") in 
slowly varying functions will be neglected. 

One also has 

PL = OW/<)J (52) 

from which 

% 

P 3 

Because only zero and first order terms are retained all other variables 
may be assumed to transform identically, except for P 2 , which appears 
explicitely in the zero-order Hamiltonian. For P 2 one finds 


(Jc/t\ s «)^ sin ( 2 fl.fl) 

(55) 

(JSe/flc)^ cos (2flfl) 

(54) 


p 2 = / dwA>q 2 = p 2 + ('aw/'os )(T>5Aacr) + o(e 2 ) 

= P 2 + AP 2 + °(£ 2 ) (55) 

It may be shown that the integral QW/'dS reduces to 


- % 


V 


2 S 


^w/7)S 


(J/4* 5 ) sin Ujtfl) 


(56) 
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wMr»h allows A p 2 to be computed. The new Hamiltonian may thus be written 
H = H (0) + H (l) + 0(S 2 ) (57) 

where the zero-order part is 

= J6J/2JT + (7<T) 2 p 2 2 /2m (58) 

and the first-order correction is 

= m(v^°^ v^ 1 ^) + P 2 V<T» (v^ + v^) 

+ p 2 A P 2 /2m (59) 

DRIFTS 

From the transformed Hamiltonian we can derive the slow variation of the 
guiding center coordinates ( U , $ ) , which may be conveniently redefined 
as (cp^/e, q 1 ) . This represents an 0( E. ) difference from ( 55 ), but since 
only the lowest order of the motion will be derived, this difference may be 
ignored. 

One can always assign to the guiding center a velocity V satisfying 
(subscript "c" omitted) 

c/ = (V*Vo() + '"do^/'Dt 

(60) 

f m (v.vj) + 

In the present case O/ot vanishes and Hamilton' s equations give 
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U - cpj/e =» -(c/e) 'dB/'i^ » -(c/e) 'dH/'^p 

p » 4 X / 7> H/d a (c/e) ^<)Bt/^«< 

Hence 


(61) 


V x B a V x ( V<*xVj3) 

a Voi. (V'Vf) - vp (V*VoO 

. a Vo^ 

a (c/e) Vj^H (62) 

The component of V parallel to B is not determined by (60) and by the 
equation that follow it and will therefore be set equal to zero, since we only 
are interested in motion perpendicular to field lines. Forming the vector 
product with B then gives 

V a (c/eB 2 ) Bx V^H (63) 

Because higher order corrections are not derived, we may s gaj n ignore the 
difference between V x and V . Using (58) we find that V H consists 
of two parts, each of which contributes to the "drift motion" represented 
by V . The first part contributes the so-called gradient drift 

Zg - (c/eB 2 ) j/2T[ (B * Vc*) ) 

a (c/eB 2 ) Ju)/2IT (B x VB) (64) 
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Because of ( 50) 

Jco/ 2 ^ * + 0(e) (65) 

from which 

V a ( v 2 /2WB ) B X VB (66) 

which is the customary expression for this drift. The second term contributes 
the so-called curvature drift 

V = (c p 2 /2emB 2 ) B X V(V<T) 2 (67) 

-curv c. — 

By (25) 

i V(v<r ) 2 . vor • vvqt 

a [1+ 0(e) ] • [vb -V«*7a -7fVb + 0( t 2 ) ] 

a £» VB + 0( £ 2 ) (68) 

with the last equality following the vanishing of B»VX and B* vji . To 
lowest order one may also replace p 2 with m v /f , where is the 
velocity component parallel to the field, giving the drift in its usual 
form as 

V = ( v 2 /«) BX(B.VB) (69) 

—curv " 

If the 'D/'T) t terms of (60) had been retained, then (63) would have 
contained an additional term 

V p » (c/eB 2 ) BX jv^ (''&<* /7>t) - V <* ( /Dt) } (70) 
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representing the velocity with which field lines appear to move due to 
'7 )o(/' 7> t and ''dJV'Dt (^)(^0(l5) # this certainly does not 

represent all effects of time dependence, a term of this form may be 
expected to appear when one deals with a time-dependent field. 

THE FIRST ORDER INVARIANT 

By (ll) and (59) the adiabatic invariant associated with gyration around 
field lines is 

J* » J + (2H/cj m ) £ia 2 (v^» v^) oac + P 2 A P 2 + 

+ m p 2 V<r . ( v ( ° J + Z (l) ) osc j + 0( £ 2 ) (71) 

where 6. is not explicitely written since it is already included in the 
present definition of and where the subscript "osc" is omitted from 

the term containing A p 2 since, by (55) and (56), this term is purely 
periodic with zero average. 

The above result is not in useful form, since J is not directly 

2 

observed. It is better to eliminate J in favor of v x , the square of 
the total velocity component perpendicular to the field, the field's 
direction being chosen as that existing at the guiding center. By (32), 

(50) and (51) 

J = 7T a ( v^) 2 / cO + 0( fc 2 ) (72) 

In 0( t- ) terms this may be used with the superscript zero omitted, as 
in (65). In the first term of (7l), however, J must be expressed accurately 


to order 6 


21 


v^ + A v x 


Then with 0( 8 ) accuracy 


= ( v ( x 0) ) 2 + 2 ( v ( ° } . A v J 

= JC0/27T + 2 ( V^. AvJ (74) 


Prom (43), (25) and (35). to lowest order 




V ( j° + m” 1 p 2 v x cr 


■ v^ - m” 1 P 2 j^ (c/e) v x a + SjVxV] (75 

Before substituting this in (71) one may note that to the order of that 
equation 

n p 2 7<T« ( v^ + v^) o np 2 (\J<r. v^) + np 2 (J. if 1 ^) 

The combination of (7l) and (75) therefore gives 

J ** (:nmv 2 /ft>) - (2JT/o)m) ^m 2 ^v^. v^^ + 

+ ■ p 2 ^yif • - ® p 2 (2 • i (1) ) 0 so “ p 2 ap 2 } 


+ 0 ( ir) 


The averaged terms vanish. In evaluating them we may discard all terns 


that are linear or cubio in components of £' 


inoluding such faotors 


as Qj , Pj and J — since sueh terms average to zero. By (75)* the tern 
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in (77) involving is linear in and vanishes upon averaging. 

In deriving the term preceding it «e nay omit all terns of 

that are quadratic in oscillating quantities* i&ich hy (44) only leaves 

» 2 < Xj?* X ( i > - • *2 < x^j >77r) (78) 

Nov it nay be shovn from (3l) that 



- - i(vf ) 2 

(79) 

1 x I 

- Vy B 2 /(Vf ) 2 

(90) 

Combining these with (32) and (39) gives 


J - (VW*>B) P 3 B 2 /&f ) 2 + 1 (e/c) Qj (V J ) 2 } 

(81) 

By (56) the average of Pj 

vanishes, vhile (compare eq. 50 ) 


(e/c) 2 <(<^y (?j» 2 - 

<P 3 2 > b 2 /(^) 2 


a 

Ju> n/2 3T 

(82) 

Therefore* from (32) and (8l) 



< y - (j/2nn 

»B) [l(V^) - (Vj*)l] 

(83) 

This is an antisymmetric dyadic and its double scalar produot with the 
symmetric W<T vanishes. Thus if one ignores correction terns purely 
periodic in A. , v^ /b is invariant even to order £. , provided 


B is evaluated at the guiding center. This result has already been 
no ted ^ 16 ^ . 
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To derive the first remaining 0( £ ) tern in (77) one requires the 
oscillating part of . From (44) 

n >2 S-2 (1 L - ^ ^ !77,r * 

+ p 2 (« - if , 2 (i- v >)) 0 sc (84) 

Actually one must excercise care with the middle term appearing on the 
right. When ? • V v ( ° } was originally derived for this tern, the V 
operator did not act on P ? and ^ hut only on the factors preceding 

them: 

m BJiVv^ » (b|jVVoOP 3 - (e/cKlE^VpCOj + >P 3 ) (85) 

(the term containing V> has a vanishing factor B* VjJ ). At the present 
stage, however, and are regarded as functions of (J,^-) 

through the relations (55) and (54). These relations also introduce a 
dependence on guiding center variables, through the function S which 
they contain. Thus when v^ is viewed as a function of ( J, a ) » 
contains extra terms with 


V Qj » - (Qj /2S ) V? 

VPj = (p 5 /2 8)VS 


( 86 ) 


In equation (84), however, this redefinition of V v ( ° } makes no diffe- 
rence, since the extra terms are there accompanied by the vanishing 
factors (| • V* ) and (| • ) • 
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By equations (25) * within higher orders, (£ • ) equals . 

Also, by (54) 

( p 3 2 ) osc ■ (e/c) (JS/2 tt) cos (4^ft) (87) 

Using preceding results and expressing A P 2 through (6) gives as the 
final form of the invariant 

J* = JTmvf/tD + (2n/<j) P 2 2 (Bf JVVr) + 

+ m p 2 (B hVv ( ^) 08C - (e/c) p 2 cos (4^a) 

+ P 2 (j/43\ % ){n>l/^(T) sin (4**1) J + 0(£. 2 ) (88) 

t 

The various first-order terms may be transformed to the forms used by 
(15) 

Northrop and by others' . This yields all the usual first-order terms, 
plus one term which has not been accounted for and which may be due to a 
ddff Qren t definition of the guiding center, due to an error or else may vanish 
by some identity. 

THE TRANSFORMED HAMILTONIAN 

The transformation which replaces ( J , £1 ) with ( J*, A.* ) affects 
the remaining variables only in ‘their second order.- If the generating 
function of; this transformation is of the form (2), then (for instance) the 
variable p 2 * replacing p 2 is (compare equation ll) 

P 2* " p 2 " 6- ^ q 2 + 0(e 2 ) (89) 

and since q 2 is "slow", the second term is also of order £ 2 . We will 
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therefore not change the notation of such variables but assume that they 
transform identically. The new Hamiltonian H* then has a zero-order 
part resembling (58) t 

H* . J*40/2JT + (V<T) 2 p 2 2 /2m + + 0( £ 2 ) (90) 

The first-order part is given by (9) and (59) as 

H*^ 1 ) » v^* (m v^+ PgVjCr)^ + P 2 + m_1 P 2 </*P 2 )> 

(91) 

but as was shown in the last section, only the middle term here does 
not vanish. The factor ( VC ) will be kept in the zero order part, even 

though it could be resolved into terms of zeroth order ( = unity ) and 
of first order. In principle, the first-order part of this expression could 
cause trouble, since it contains terms proportional to oC and , and 
one of these may be angle-like and grow without bounds. Formally one 
can resolve this by switching to a new ( «< , ^ ) system when the 
angle-like variable has completed a full circuit; in practice this ambi- 
guity may be expected to disappear when results are expressed in observable 
quantities and in any case, at a later stage (V0") will be transformed 
away. 

Using (44) one finds, to the order of the approximation 




£ m"*S> 2 '^V'* <r 

(92) 
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where the contribution of the middle term is omitted since it is linear 
in oscillating quantities. 

(l7) 

The evaluation of the first term requires the vector identity 

A B : V C = BAsVC + (Bxa).(V*C) ( 93) 

(the notation of the preceding equation and of the two that follow is 
independent of the one used in the rest of this work). How let A and 
C be two orthogonal vectors. Then 

(A.C) = 0 = A* VC + C* VA + Ax(Vxc) + C*(V*A) (94) 

Performing the scalar multiplication of (94) with an arbitrary vector B 
and applying (93) then gives 

A B : VC = - C B : VA (95) 

By (39) the magnetic field and 8 are orthogonal; therefore, using (83) 
and (93) 

<^B ?: V x v^)> => - :7B 

= - ( J* / 2JT mB ) [l - 7f l] : 7B 

= - ( J* / 2JTmB ) (vj& x T) • (Vx B) (96) 

Here one must use J as variable, since we are dealing with a term of 
H* . By (31) 

Tx7f = 


B 


(97) 
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so 

(l? ! V2^) - ( J*/ 2Jrm ) B • (v x B) (98) 

To lowest order J* equals J , and by (82") one may therefore 
substitute in the last term of (92) 

<P 5 2 > - J*Se/2Jrc + 0(£) (99) 

giving as final result 

H*^ = (p 2 J*/21TmB) B • (VxB) - (p 2 J*$ e/4JTmc) 't>A/a<T (lOO) 

The possible applications of H* include the associated Liouville 
equation, which should lead to a gyration-averaged form of the kinetic 
equation for a collisionless plasma, of the type first derived by Chew, 
Goldberger and Low^^. 

THE SECOND PERIODICITY — THE LIMIT £ = 0 

Consider the limiting case of (90) in which higher order terms and 
the dependence of on (p^, q^) may both be neglected. The Hamiltonian 

H* (0) = (q 2 )/2 + (VO-) 2 p 2 2 /2m (lOl) 

may then be viewed as describing the motion of a particle in a potential 
proportional to a>(q 2 ) . If the dependence of cd on q 2 (i.e. on 
of the guiding center) has the form of a potential well deep enough to 
trap the particle, the motion will be periodic and one can again derive 
action-angle variables and an adiabatic invariant. 

The definition of <T is not unique, since there exists a large 
arbitrary element in the choice of ( o( , 0 , s ) , but the dependence 
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of field quantities on it always describes the variation along a field line, 
since the other guiding-center ^variables already specify the field line 
itself. If the field intensity B (and therefore W ) is large at two 
separated points on a field line and lesser between than, a "potential well 
geometry" will exist and under suitable conditions the particle will be 
trapped in it. These conditions involve the energy integral E of (lOl) : 
since the second term in that equation is non-negative, if the particle 
starts between two values of q 2 for which J*o) (q 2 )/ 2 JT exceeds E , it 
will be trapped. 

Assuming that periodicity exists, we now seek the generating function 

¥ ( j 2 » ^ 2 } to new variables (^ 2 *^ 2 } that are action-angle variables of 
(lOl). Since this is the limit £ = 0 , both p^ and q^^ are to be con- 
sidered as constant parameters that do not vary in time; independently of 
£. , J is always a constant at this stage and may be removed from consi- 
deration as a canonical variable. 

If J*) is the period of this motion, XL 2 must increase 

by unity during each period and therefore the new Hamiltonian must depend 
on J 2 through a term J 2 / T . To derive T one must solve the zero- 
order motion in the given potential well, something that can often be only 
done numerically or approximately. This is one point of difference between 
this periodic motion and the gyration around field lines discussed earlier; 
with the gyration the basic frequency requires only knowledge of the local 
magnetic field, while here a more complicated solution of the zero-order 
periodicity is required. 
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Assuming T to be known, the Hamilton-Jacobi equation becomes 
J*w/2J1 + (v(T) 2 (^W/^q 2 ) 2 /2m « J 2 / T(p r q r J*) + F (102) 


where P is a yet undetermined function of the slow variables (not related 
to eq. 26 ). Prom this 


(2m)‘ 


[j^T + F - J* 0 )/ 2 JT] * [dq g Av<r) 2 ] (103) 


giving 


(V(T) 2 p 2 2 /2m = Jg/T + P - J*6,)/2JT 


(104) 


which confirms the form of the new Hamiltonian as the one given on the 
right of (102). Also 


n 


2 




(•/*)* 



+ P 


(105) 

J*U)/2 JT] " * [ dq 2 /( V<T) 2 J 


The significance of this is easily seen. Since J*6)/2JT is the energy 
associated with the component of the motion orthogonal to the field, we 
have (at least to zero order) 


= P 5 / n 


Also, since along a field line ( ) do not change 


dq 2 = 


d<r = ds + 0( g, ds ) 


(106) 
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Hence, if X is the time of motion evaluated using v y/ alone 


n 


T 


,-i 


(ds/v.) = X/ T 


( 107 ) 


In order to express T we use the fact that CL . increases by unity 

(105) 2 

each period. Integra ting^over one full period then gives 


r 


T = m O (dq^p^ 




(ds/y,,) 


(108) 


It is interesting to note that the familiar formula 


f 


'2 - J p 2 dq 2 (109) 

is not obtained a3 part of the calculation. It enters as follows 
since Cl increases by unity each period 


f 


d n 


^ dq 2 
<j> OV^^) dq 2 

(^) /7) J 2 ) j> (()¥/Dq 2 ) dq 2 


(7) / 7 > Jg) p 2 dq^ 


( 110 ) 


The last condition is satisfied if (109) holds, but more generally it 


only requires J 2 to have the foxm 


^ P 2 d 9 2 + K(q x , p x ) 


( 111 ) 
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The function K is arbitrary and varying the choice of this function 
(with H remaining constant) results in different choices of the function 
F in (102) and in the equations that follow there,, 


However, if we wish to use the transformation (l03) as the starting 
point for a perturbation expansion there exists an additional considera- 
tion which die tates t ^ e choice of J 2 satisfying 109) : unless such a choice 
is made, the transformation of the subsidiary variables (p^, q^) for 
finite £ will contain a secular term. This point will be discussed 
in the next section; the function F to be used is than implicitely 
defined by the relation 


J 


2 



J*u)/2Jt] * [dq^CvcT) 2 ] 

( 112 ) 


TOB SECOND PERIODICITY — FINITE £ 

To transform (90) with (p^, q^ (but not J*) restored to canonical 
status, equation (103) must be modified to 

W(J 2 > ^ 2 * ^ 2 / 9^) = * 0 + P 1 ^1 

where (P^, Q^) are new variables (though older notation is used in 
labeling them) and Wq is the expression on the right side of (103), 
evaluated at (cP^/e , q^ » * 

The transformation of (p 2 » q g ) resembles the one previously 
obtained — in particular, (104) is again obtained. Because T, F and u> 
are now evaluated in "mixed" variables, this confirms the form of the 


zero-order Hamiltonian 
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H = J 2 /T(p r V J*) + p ( p i» V J *) + o(e) (114) 

only if the transformation of the subsidiary variables (or at least of q^) 
is a near-identical one. For q^ we have, by (113) 

Qj. - 

The second term here contains which lowers the term by one 

order in 6 — P^ being slow — but it also contains Wq which is an 
open-ended integral in q 2 . Unless this term is either periodic or 
independent of P^ , it may grow without limit and after 0( €. -1 ) periods 
this growth may be large enough to offset 'd/'O P^ and give a transfor- 
mation which no longer is of near-identity type. 


Let q 2 (k) denote the value of q 2 at the end of the k-th period 
and let N be a large integer, of order . Suppose the system is in 

it3 (N + l) period: then 



The last term in (115) is 0(l) and differentiation reduces it to 0(£) , 
allowing it to be ignored. The summation term on the other hand is large: 
it consists, of two parts, the first of which does not survive differenti- 
ation, since it is equal to NJ 2 . The second part could raise (115) to 
0(l), but if J 2 is defined as in (109) K vanishes and this part does 


not arise. 
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It may be remarked here that the 3ame consideration is implicit in 
the earlier transformation (Py Q^) — ^ ( J, H ) , generated by (46). In 
that case, again, expressions such as 'DW/'Dp^ appearing in the 
transformation of subsidiary variables could give rise to secularly 
increasing terms if the open-ended integral ’"DW/'DS exhibited such 
growth. As is shown in (56) this integral is in fact periodic in A and 
such growth is avoided. 


The first-order Hamiltonian is 


H = J 2 /T + F + + 0(& 2 ) (H6) 

where is obtained from (lOO) by replacing (p 1# q^ with (P^ Q x ) 

and eliminating (p 2 » q 2 ) *>y means of (104) and (105). The replacement 

of q 2 by n 2 is rather difficult, but fortunately it may be by-passed 
in evaluating the first-order adiabatic invariant. By (ll) and (8), this 
invariant is ( E not written, averaging taken over f\~ 2 ) 


J 2 + 


T ( - <H UJ > ) + 0(a 2 ) (117) 


,(D 


Let q 2 be retained in as an auxiliary variable replacing jfL 2 . 
To find <( one introduces q 2 as integration variable instead 
of Jfl 2 , using (105) (other variables may be ignored to lowest order): 


<H U) > 


I H (l) dfL; 


s 


H (l) (DA/t^) dq 2 


r (118) 
(m/2)^ T" 1 <P [j^T + F - J # £0/2 7r] “ * [iq^?^) 2 ] 
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No attempt has been made to compare (117) and (118) to the first-order 
derivation of J* presented elsewhere^ 2< ^ . 

THE THIRD PERIODICITY 


Assume now that the transformation eliminating D- 2 has been carried 
out, so that the new Hamiltonian contains J* as a (constant) canonical 
variable describing the longitudinal motion. In this transformation 
and Q 1 also undergo near-identity change and their transformed versions 
will be simply denoted P and Q , subscripts no longer being necessary. 

The new Hamiltonian is 

H(P, Q, J*. J*) » J*/T(P, Q, J*) + p + <H (1) > + 0 (e 2 ) 

(119) 

where <^H^ ^ and P are the same as in (116) and (118) but with 
new variables replacing old ones, a substitution which affects only higher 
orders since this again is a near-identity transformation. Again this may 
be used, via Liouville's theorem, to derive a kinetic equation for a collision- 
less plasma trapped in a mirror-like geometry, with averaging over both 
gyration and longitudinal motion. It may be mentioned here that the 
"Hamiltonian properties" of J* / T were already noted by Northrop 

The Hamilton-Jacobi equation of (119) now seeks the generating function 
W(p» Q) of a canonical transformation to a new set (p, q) in which p is 
constant and its conjugate q is linear in time. This function satisfies 


H [(WS«), 0, J*, J*] 


Vp 


( 120 ) 
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with V a yet undetermined, constant which could involve J* and J* . 
Assuming that 'TlW/'&Q can be extracted from this 

^W/aQ = f(Vp, Q, J*, J*) (121) 

the equation is easily solved by quadrature. 

The physical significance of this somewhat formal procedure hinges on 
the fact that independently of the foxm of W , a canonical transformation 
generated by it satisfies 

'Mp, q)/D(P, Q) = 1 (122) 

This is another way of saying that a canonical transformation represents 
an area-preserving mapping in the (p,q) phase plane, a fact already noted 

(t) 

by Gardner' ' and by many others. One may prove (122) by direct substitution. 

Now P and Q themselves are related by near-identity transformations 
to the variables and of (35) which, apart from an unimportant 
factor (c/e), carry the connotation of guiding center Euler potentials. 

These potentials, however, are far from unique: if ( o( , ^ ) represents 

, - i 

one choice for them, equivalent choices oC ( ol , p ) and y ( d , JS ) 
may be used to describe the same magnetic field, provided 

°C,P ) = 1 (125) 

In view of (122) the new variables (p, q) may be identified as 
describing ( to lowest order) a different choice of ( , ^ ) for the 

field in which the particle moves. One may say that among the many different 
choices of { ^ , f ) describing the given field, there exists one choice 
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— or rather, a set of choices — such that a given particle with speci- 
fied J* and J* » on the average (average over XL and X\ 2 )_, stays 
on the same ^ and advances uniformly in . The Hamilton-Jacobi 
equation (120) then provides a method for deriving Euler potentials 
belonging to this particular set. 

If the new ^ — that is, q — is an angle-like variable, the motion 

will again be periodic. This will always occur in axisymmetrical fields 

(if the other two periodicities exist) for if then is initially chosen 

of 

as the azimuth angle Y , it will be absent from the Hamiltonian in all/iits 
transformations. Thus in this case Q does not appear in (119) and the 
last transformation is trivial. However, the field needs not be axisyrnsnet- 
rical for the third periodicity to occur, as is demonstrated by the motion 
of radiation belt particles in the earth's magnetic field, which may be 
viewed as a distorted dipole field. 

If the motion is periodic it will possess an action variable 

J 3 = P dq (124) 

(22) 

This is commonly known as the third invariant or the flux invariant , 
since if p and q are viewed as averaged Euler potentials, then 

= (e/c) ^ (125) 

and the integral appearing here equals the magnetic flux ^ embraced by 
the particle's orbit in one period. From general theory one ^ 
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expect it to be adiabatically conserved under slow perturbations. Our 

field cannot provide such perturbations, since all variables have 

( 22 ) 

already been accounted for; it may be shown, however , that will 
be adiabatically conserved if the field varies slowly in time or if an 
electric field is slowly applied to ito 
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APPENDIX 

We wish to relate the first-order correction of (88 ) to geometrical 
quantities. Extracting a factor B from the first tern and using ( 93 ) 
gives 

1 J* VVs Q = $ B s VVs Q 

= IS i ( VB - 7<^7a - Vf 7b) + 0(8 2 ) 

= |B:VB + 0 (£ 2 ) (a- 1 ) 

Forming the dot product of $ with 
V*(BxB) = 0 

= B*VB - B- VB + B ( V • B ) (A-2) 

and using ( 93 ) once more gives, finally, 

B$ J V V s Q - 8 B s V B 

=» B £ : V B - (v x B)/cO (A-3) 

which, apart from non-geometrical factors, equals the v ^ 2 term given 
by Northrop ^ 

Before evaluating the re ma i n ing terms, all of which are proportional to 
Pp to ) » it is useful to derive V * v^x^ » i * 1 the system of 

variables that includes (j, fl ). We have, by ( 31 ) and (86 ) 
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a (V x = V * I p 5 Vo< - (e/cXQ^ + > P^) vf>] 

= ( V P 3 x VoC ) - (e/cXVQj + W P 3 + P 3 V> ) * ^ 

. (P /2S )(VS x T) + (e/cX^M KvS^vf) - (e/c)P 3 (V^7f) 

(A-4) 

Using (32) and (97 ) 

(m 2 /B) v ( °^ (V X £ 0) ) = (m/B) T - (e/c^Vf]* ( V x v^) 

= - (% p 3 A )(B * VS) + (e/c) P 3 2 (B*V>) 

(A-5) 

By (41) 

(B'VS) = QS/^cr + 0(E 2 ) (A-6) 

Comparing this to (84), (77) and (56) shows that the 0( £ ) correction term 

. ,* . 

proportional to v f/ in J xs 

(2K/o)) P 2 { 1 ^ - (1/2W) v ( ° } -(v x v ( ° } )} Qac U-7) 

Replacing the factor by w(i> B) and applying ( 93 ) gives for 

this term . . 

(JT/oj) P 2 {iP : V v ( x 0) + JB :V v ( x ° j QSC (A-6) 

The average of the first term in (A-S) is given by (98 ) and is easily 
subtracted; the result agrees with Northrop' s v //f term. The second term 

remains unaccounted for and has not been further evaluated. 
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